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ABSTRACT 



Context. The most nearby active galaxy Cen A has attracted considerable attention as a detected TeV gamma-ray and possible ultra- 
high energy (UHE) cosmic-ray emitter. 

Aims. We investigate the efficiency of particle acceleration close to the supermassive black hole (BH) horizon assuming that accretion 
in the innermost part of the disk occurs in an advection-dominated (ADAF) mode. 

Methods. We analyze the constraints on the achievable particle energies imposed by radiative losses and corotation for conditions 
inferred from observations. 

Results. We show that for an underluminous source such as Cen A, centrifugally accelerated electrons may reach Lorentz factors of 
up to y ~ (10 7 - 10 8 ), allowing inverse Compton (Thomson) upscattering of ADAF sub-mm disk photons into the TeV regime with 
an associated maximum (isotropic) luminosity of the order of a few times 10 35 erg/s. Upscattering of Comptonized disk photons is 
expected to lead to a TeV spectrum L v oc v~" c with a spectral index a c - (1.5— 1.9), consistent with H.E.S.S. results. The corresponding 
minimum variability timescale could be as low as r^/c ~ 1 hr for a typical light cylinder radius of - 5 r s . While efficient electron 
acceleration appears to be well possible, protons are unlikely to be accelerated into the extreme UHECR regime close to the central 
black hole. We argue that if Cen A is indeed an extreme UHECR emitting source, then shear acceleration along the kpc-scale jet could 
represent one of the most promising mechanisms capable of pushing protons up to energies beyond 50 EeV 
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1. Introduction 

Centaurus A (NGC 5128, Cen A) is the closest (distance d ~ 3.4 
Mpc; Israel 1998) and one of the best studied active galaxies. 
Optically, Cen A is an elliptical galaxy undergoing late stages 
of a merger event with a small spiral galaxy. Radio observa- 
tions detected a complex FR I morphology with a subparsec- 
scale jet and counter-jet, a one-sided kiloparsec-scale jet, two 
radio lobes and extended diffusive emission. VLBI observations 
of the subparsec-scale jet indicate that Cen A is a non-blazar 
source with a jet inclination angle i ~ 50° (Tingay et al. 1998), 
or perhaps somewhat smaller if large-scale observations are reli- 
able tracers (Hardcastle et al. 2003). The center of its activity is 
a supermassive black hole (BH) with a mass inferred to be in the 
range ot B h - (0.5 - 1.2) x 10 8 M G (Marconi et al. 2006; Haring- 
Neumayer et al. 2006), corresponding to a Schwarzschild scale 
r s a (1.5 -3.6)x 10 13 cm. 

Cen A is the only AGN of the non-blazar type detected at MeV 
(COMPTEL: Steinle et al. 1998) and GeV energies (EGRET: 
Sreekumar et al. 1999). The nuclear spectral energy distribution 
(SED) of Cen A, as inferred from non-simultaneous data, ap- 
pears to consist of two peaks, one reaching its maximum at sev- 
eral times 10 13 Hz and one peaking around 0.1 MeV (Chiaberge 
et al. 2001; Meisenheimer et al. 2007). The SED below 1 GeV 
has been successfully modeled within a simple jet synchrotron 
self-Compton (SSC) framework, assuming Cen A to be a mis- 
aligned BL Lac object (Chiaberge et al. 2001). In the early days 
of gamma-ray astronomy, a tentative detection (> 4<x) of Cen A 
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at very high energies > 0.3 TeV was reported during a phase of 
high X-ray activity in 1972-1974 (Grindlay 1975). H.E.S.S. re- 
sults (Aharonian et al. 2009a) have established Cen A as a TeV 
emitting source (the second radio galaxy after M87) with inte- 
gral VHE flux of ~ 1 % of the Crab Nebula. 
Chandra and XMM-Newton observations between 2 and 7 keV 
indicate that the nuclear X-ray continuum spectrum in Cen A 
may consist of both a disk and a jet component, the disk contri- 
bution being consistent with a hybrid disk configuration where a 
standard disk is truncated at r, and replaced by an ADAF in the 
inner regions close to the central black hole (Evans et al. 2004; 
cf. also Pellegrini et al. 2005 and Meisenheimer et al. 2007). A 
transition to an inner ADAF disk could also explain the lack of 
a big blue bump UV feature expected in the standard disk sce- 
nario (Marconi et al. 2001). The accretion rate for Cen A seems 
uncertain by about an order of magnitude, with model-dependent 
estimates ranging from several 10~ 4 to some 10~ 3 mEdd- Taking 
the measured nuclear X-ray luminosity of L x ~ 5 x 10 41 erg/s 
(Evans et al. 2004) as a robust upper limit to a possible ADAF 
disk contribution gives m ~ 0.004»iEdd, which for the canonical 
viscosity parameter a ^ 0.3 is still below the critical mass ac- 
cretion rate m cr j t - 0.3 a 2 m£dd required for a two-temperature 
ADAF to occur (Narayan et al. 1998; Yi 1999). 
If the first SED peak below 10 14 Hz is indeed caused by 
synchrotron emission from the jet (Chiaberge et al. 2001; 
Meisenheimer et al. 2007), inverse Compton upscattering in a 
conventional homogeneous SSC framework is unable to account 
for VHE y-rays in the TeV regime. On the other hand, it could 
be that the situation for Cen A is similar to that for M 87, where 
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efficient particle acceleration close to the BH horizon seems to 
be responsible for the production of the observed TeV y-rays 
(Aharonian et al. 2009b; Rieger & Aharonian 2008a, henceforth: 
RA08a). Here we analyze possible VHE characteristics associ- 
ated with this scenario. We show that for parameters inferred 
from observations, inverse Compton upscattering of ADAF disk 
photons by centrifugally accelerated electrons could result in 
variable (~ 1 hr) VHE y-rays with a relatively hard TeV spec- 
trum. We show that the mechanism considered is unable to ac- 
celerate protons to ultra high energy cosmic-ray (UHECR) ener- 
gies, although shear acceleration along the jet could potentially 
provide a means. 

2. Cen A as a TeV y-ray source 

2.1. Particle acceleration close to the BH horizon 

Magnetic fields are widely considered to play a key role in the 
formation of collimated relativistic outflows from compact ob- 
jects. According to, e.g., MHD scenarios, magnetic flux dragged 
in by the accretion flow and suitably amplified by dynamo ac- 
tions in the inner accretion disk can build up a rigidly rotating 
(quasi-dipolar) disk magnetosphere. Depending on the field-line 
rotation and the current distribution, the bulk of the plasma may 
then be centrifugally accelerated along open flux surfaces to rela- 
tivistic speeds (F/, ~ 10) and become collimated outside the light 
cylinder r L (5 - 10) r s (e.g., Camenzind 1995; Fendt 1997). 
In what follows, we explore an approach where efficient parti- 
cle acceleration is supposed to occur close to the light cylinder 
(e.g., Machabeli & Rogava 1994; Gangadhara & Lesch 1997; 
Rieger & Mannheim 2000; Osmanov et al. 2007). Within the 
MHD framework this could occur if, e.g., the longitudinal cur- 
rent becomes smaller than a critical value (Beskin & Rafikov 
2000). Let us suppose, thus, that a significant fraction of the 
electromagnetic energy may be transformed into the kinetic en- 
ergy of particles by means of magneto-centrifugal effects close 
to tl. Consider for illustration a charged test particle, corotating 
with the magnetic field and gaining energy while moving out- 
wards (Machabeli & Rogava 1994). The motion of this particle 
may be analyzed conveniently in the framework of Hamiltonian 
dynamics: For length scales much smaller than the curvature ra- 
dius of the field, the field topology is monopole-like and the as- 
sociated Hamiltonian H for a particle of rest mass mo moving 
along a relativistic ally, rigidly rotating field-line (angular veloc- 
ity Q = c/rL = constant) becomes a constant of motion 

H = ymoc 2 (l - r 2 /r 2 ) — const. , (1) 

where y = (1 — r 2 /r 2 - r 2 /c 2 )~^ 2 is the particle Lorentz factor. 
Thus, as a particle approaches rL, its Lorentz factor might in- 
crease dramatically (provided corotation holds; however, it can- 
not become infinite, see below). Using Eq. ([TJ and the definition 
of y, this can be expressed in terms of a characteristic (minimum) 
acceleration timescale (cf. RA08a) 

_ y 1 

?acc ~ t * 2Qm!/4 r i/2 ' (2) 

valid for y » 1, where m = l/(y? [1 - r?/r 2 ] 2 ) depends on the 
initial injection conditions. 

2.2. Constraints on achievable Lorentz factors 

Under realistic astrophysical conditions, the acceleration ef- 
ficiency is constrained by radiation reaction (e.g., inverse 



Compton in the ambient photon field), the breakdown of the 
bead-on-the-wire (BW) approximation for a single particle or 
the bending of the field line with increasing inertia. For a mass 
accretion rate in the concordant range (2 - 20) x 10~ 4 mEdd 
(e.g., Marconi et al. 2001; Evans et al. 2004; Pellegrini 2005; 
Meisenheimer et al. 2007), the ADAF equipartition magnetic 
field fi eq °c m B Jj /2 r- 5 4 m 1/2 (Yi 1999) becomes B eq (3 r s ) ~ 500 
G, suggesting (radial) jet magnetic field strengths near the light 
cylinder of B{ri) ~ 50 G. In the case of BL Lac objects, transi- 
tion radii from tens to several hundred r s have been inferred from 
observations (e.g., Cao 2003; Wu & Cao 2005). Similarly, the 
apparent lack of a standard disk feature in the observed SED of 
Cen A seems to indicate as well that r, > 10 r s . By assuming that 
all of the measured X-ray luminosity in Cen A of L x ~ 5 x 10 41 
erg/s is produced by an inner ADAF, the implied energy densi- 
ties is f/ph ~ L x /(4nrfc) ~ 1 erg/cm 3 for r, ~ 30 r s . A similar 
constraint might be obtained based on the ADAF synchrotron 
peak luminosity Lr ~ v s L Vs oc m^m^v 7 / (Yi 1999) close to 
rL, where v^l) - 7 x 10 11 Hz, which give Lr ~ 3 x 10 39 erg/s, 
using the fiducial parameters m = 5 x 10~ 4 mEdd, '«bh = 10 8 M o , 
and ^ ~ 5 r s (giving B(ri) ~ 100 G), implying a photon energy 
density t/ p h ~ 0.4 erg/cm 3 . 

(i) Inverse Compton losses (Thomson regime) for electrons, oc- 
curring on timescales of f 000 i 3 x 10 7 /[y i/ p h] s, then imply (see 
Eq. d2]i) that achievable electron Lorentz factors are limited to 

^^^(^(i^) 2 . (3) 

(ii) For motion along the field synchrotron losses are negligible, 
while curvature losses, occurring on characteristic timescales of 
t c 180 R 2 ,(mo/m e ) y~ 3 , do not impose stronger constraints 
(compared to Eq. (O) for curvature radii R c ~ r^. 

(iii) The validity of the BW approximation for protons, on 
the other hand, requires that their characteristic acceleration 
timescale is always longer than the inverse of the relativistic 
gyro-frequency (cf. Rieger & Aharonian 2008b), which con- 
strains the achievable proton Lorentz factors to be 

/:p ,^^(^fm 2/5 . (4 ) 

/max ' p \50G) \10 14 cm/ v ' 

2.3. Expected TeV characteristics 

The above considerations suggest that electron acceleration 
close to the central BH in Cen A may produce Lorentz fac- 
tors y ~ 10 7 . This would allow Thomson inverse Compton 
upscattering of ambient disk photons to TeV energies of e ~ 
ym e c 2 - 5 (y/10 7 ) TeV. The resultant TeV spectral index at the 
high energy end then traces, over a decade or so, the underlying 
disk seed photon spectrum (Rieger & Aharonian 2008b). In the 
ADAF case, comptonization of thermal synchrotron soft photons 
is expected to lead to a power-law-like tail in the disk spectrum 
above v s , i.e., L v =s L v ^(v/v s y ac , where a c = -lnT esc /lnA de- 
pends on the accretion rate (Mahadevan 1997; Yi 1999). For the 
inferred mass accretion rate of Cen A, m (2 - 20) x 10~ 4 WEdd, 
the spectral index becomes a c (1 .2 - 1 .9) assuming a standard 
ADAF viscosity coefficient a = 0.3. This suggests that the up- 
scattered TeV spectrum in Cen A might in fact be softer than that 
observed in M87 (see Aharonian et al. 2006). A spectral fit to the 
VHE (> 300 GeV) data obtained by H.E.S.S. indeed suggests a 
power-law spectral index in the range 1.7 + 0.5 (Aharonian et al. 
2009a). A soft spectrum with, e.g., a c = 1.8, may then indicate 
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that accretion occurs at the lower end of the concordant inter- 
val noted above. Since particle acceleration occurs close to the 
light surface, the expected minimum TeV variability timescale is 
~ tl/c ~ 1 hour. 

Is it possible that TeV photons, produced by inverse Compton 
upscattering, can escape from the vicinity of the central black 
hole in Cen A? In principle, photons of energy E [TeV] will in- 
teract most efficiently with target photons in the infrared regime, 
i.e., of energy eiR ^ (1 TeV/E) eV. For the radio luminosity Lr 
estimated above, the implied infrared (1 eV) ADAF disk lumi- 
nosity is on the order of Lir ~ 3 x 10 39 (300)'^ erg/s. The 
optical depth for yy-absorption T yy ^ (LiR<x yy ) / ' (AjiRi^€i^c) thus 
becomes 

TrrWM ~ 3o^T ( 3 x iXrg/s) fe) (l4v) ' (5) 

indicating that TeV photons may well escape from the vicinity 
of the black hole provided the spectrum is sufficiently soft (i.e., 
for a~ 1 .4). The nuclear (pc-scale) region of Cen A, on the other 
hand, is known to be relatively mid-infrared-bright (with lumi- 
nosity of ~ 7 x 10 41 erg/s). Yet, if this flux is indeed dominated 
by emission from a dusty torus of size J 0.1 pc (see Radomski et 
al. 2008), VHE gamma-rays would also be able to escape unab- 
solved from the nuclear region of Cen A. 

One can estimate the maximum possible TeV luminosity out- 
put Lie ~ P n P c AV based on the Alfven corotation condition 
B 1 l%n > n~ym e c 2 , which provides an upper limit for the acceler- 
ated particle number density close to r^. By noting that for the 
relevant volume one has AV = 4-7r{yQ/y)r^, that P c can be ap- 
proximated by the single particle Thomson power, and that par- 
ticles are accumulated for p = t coo \/t acc > 1 (correcting a typo 
in RA08a, eq. 12), the maximum possible inverse Compton TeV 
luminosity is of order 




which is consistent with the H.E.S.S. results of Aharonian et 
al. (2009a). Since Lie increases almost linearly with black hole 
mass, a source needs to be sufficiently massive to become de- 
tectable by current ground-based instruments. In fact, under 
normal conditions (i.e., for moderate injection Lorentz factors 
yo ~ 10), the inferred TeV output for Cen A is close to the sensi- 
tivity limit (1% of the Crab flux, with 5cr in 50h) of the H.E.S.S. 
array. Variability studies on hour timescale will thus only be pos- 
sible with the next-generation CTA telescope project. 

3. On Cen A as a possible UHECR proton source 

According to the Pierre Auger results, 4 out of 27 events de- 
tected above 57 EeV may be associated with the location of 
Cen A (Abraham et al. 2007, 2008). If some of these cosmic ray 
events were indeed associated with Cen A, it seems unlikely that 
their energization could occur close to the central BH. Neither 
a Blandford-Znajek-type nor a centrifugal disk wind scenario 
seems to allow for that: For example, if the black hole in Cen 
A rotates and is embedded in a magnetic field of equipartition 

strength B Q ~ 10 4 (m/10- 2 m Edd ) 1/2 (10 8 M o /m BH ) 1/2 G, it could 
act as a unipolar inductor (membrane paradigm) and induce a 
voltage difference of the order of 

\10 8 M o / \10 4 G/ v ; 



However in reality, it seems unlikely that protons might be able 
to fully tap this potential. Even if the ordered magnetic field 
component could be as high as Bo, screening of the electric field 
in the plasma-rich AGN environment and curvature losses, for 
example, are readily expected to reduce achievable proton en- 
ergies to below ~ 10 19 eV (cf. Levinson 2000; Aharonian et 
al. 2002). Moreover, as an FR I source, Cen A is most likely 
characterized by a Kerr parameter a well below one (e.g., Daly 
2009). Centrifugal acceleration of particles in the rotating disk 
magnetosphere, on the other hand, is limited by the condition 
for single particle corotation, cf. Eq. ©. While this works well 
for electrons, efficient centrifugal acceleration of protons to en- 
ergies far beyond 10 17 eV is unlikely to occur in the vicinity of 
the central BH. We finally note, that the required Lorentz factors 
of y p ~ 5 x 10 10 (if assumed to be protons) would already imply 
a gyro-radius r gyro » 1.5 x 10 15 (y p /5 x 10 10 ) (100 G/S(r L )) cm, 
exceeding the characteristic light cylinder scale for Cen A. 
Furthermore, if radio observations are a reliable tracer of the true 
fluid speeds in the jet of Cen A, internal shocks are also probably 
not able to account for UHECRs beyond 50 EeV. On both VLBA 
subparsec-scales and VLA scales of hundreds of parsecs, the jet 
bulk speed seems to be at most mildly relativistic, not exceed- 
ing 0.5c by much (Tingay et al. 2001; Hardcastle et al. 2003). 
If the shocks are internal with characteristic speeds compara- 
ble to the relative velocity between layers, shock speeds B s well 
below 0.5c are expected. Balancing the minimum timescale for 
non-relativistic shock acceleration (Bohm limit, e.g. see Rieger 
et al. 2007) with the timescale for cross-field diffusion then im- 
plies, that possible cosmic ray energies are limited to £ m „ ^ 
Ze B r,B s , where r, is the (transverse) source size and B the lo- 
cal magnetic field strength. Even if one would allow for some lo- 
cal magnetic field amplification by a factor p (e.g., due to shock 
compression or shear-driven instabilities) and assume that radia- 
tive losses could be neglected, possible proton energies are un- 
likely to exceed 

^<-2xl0'V.(Jy (J^)eV. ,8, 

where the strength of the local field at r ~ r,/aj (where aj is 
the jet opening angle) has been expressed as B ^ pBo (rjr), 
Pa - p/4-, with e.g., p4 = 1 for a strong shock-compressed down- 
stream region. While it seems thus unlikely that protons can be 
accelerated to 50 EeV and beyond (at least if no dramatic flaring 
activity is assumed, cf. Dermer et al. 2009), acceleration to the 
ankle at ~ 4 x 10 18 eV might well be possible. 
Hardcastle et al. (2009) recently proposed that efficient stochas- 
tic particle acceleration could occur in the giant radio lobes of 
Cen A on scales of R ~ 100 kpc. In the Bohm limit, the charac- 
teristic acceleration timescale for stochastic second order Fermi 
acceleration is of order f acc ^ (c/vA) 2 (r gyro /c) (e.g., Rieger et al. 
2007). Comparing this timescale with that for cross-field diffu- 
sion, i.e., fdf ~ 3R 2 /(cr gyro ), suggests that proton energies are 
limited to 

£ raax -1.6xl0 19 p^) {—?—) {— ^-)eV. (9) 
\0.1 cl \ 100 kpc/ U0- 6 G/ 

Hence, if relativistic turbulent wave speeds can be achieved (of 
order va — 0.3c, as assumed in Hardcastle et al. 2009), proton 
acceleration to energies 5 x 10 19 eV could become possible. Yet, 
whether these conditions exist, seems unclear (cf. O' Sullivan et 
al. 2009). If (part of) the observed X-ray emission from the giant 
lobes is produced by thermal plasma emission, the implied ther- 
mal plasma density would be of order n± ~ (10~ 5 - 10~ 4 ) cm -3 , 
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i.e., much higher than conducive to fast Alfvenic turbulence. For 
typical equipartition field strengths of B ~ 10~ 6 G, the character- 
istic Alfven speed va — Bj ^An n t h m p ~ c/500 would be much 
lower than required. 

Based on the above considerations, it seems questionable 
whether protons could indeed be accelerated to energies be- 
yond 10 19 eV in Cen A. However, one possible mechanism that 
may help particles to reach beyond this limit could be shear 
acceleration along the large-scale jet (Rieger & Duffy 2004). 
Strong internal jet stratification (e.g., a fast spine surrounded 
by slower moving layers) is naturally expected for jets inter- 
acting with their environments or launched by different pro- 
cesses (e.g., disk and black-hole driven wind), and indeed seen 
in hydrodynamical simulations (e.g., Meliani & Keppens 2007; 
Aloy & Mimica 2008). In the case of Cen A, there are some 
observational hints, e.g., limb-brightening in the X-ray jet, the 
higher energy emission being more confined to the jet axis, 
longitudinal magnetic field polarization in the large-scale jet 
(Kataoka et al. 2006; Kraft et al. 2002; Hardcastle et al. 2003) 
that might be indicative of internal jet stratification. Energetic 
particles scattered across such a shear flow can sample the ki- 
netic difference in the flow and be accelerated to higher en- 
ergies. The timescale for particle acceleration in a continuous 
shear flow is inversely proportional to the particle mean free 
path f acc oc l/A (Rieger & Duffy 2004). In order to compete 
with radiative losses, and in the case of Cen A with advection 
along the jet, high energy seed particles are thus needed to en- 
sure efficient particle acceleration. Assuming a gyro-dependent 
particle mean free path A ~ r gym , the minimum acceleration time 
scale in a linearly decreasing velocity shear profile is on the or- 
der of f acc a 3(Ar) 2 c/[A y 4 (Av z ) 2 ] ~ 6 (Ar) 2 /(r gyro c) assuming 
v, ~ 0.5c (see Rieger & Duffy 2004). Advection along the jet 
occurs on a characteristic timescale t a ~ lj/v z ~ 10 12 sec, where 
lj ~ 4.5 kpc is the (projected) X-ray jet length (Hardcastle et 
al. 2007). The observed kpc radio jet slightly expands and has a 
transverse width of ~ 1 kpc close to its end (Kraft et al. 2002). 
Protons injected at the crossover to the kpc-scales, for example, 
would thus need seed energies of E ~ (1 - 5) x 10 18 eV to en- 
sure that shear acceleration can operate efficiently (using, e.g., 
Ar ~ 0.2 kpc, B ~ 10~ 4 G). Such energetic seed particles may 
well be provided by e.g. shock acceleration. In this case, protons 
escaping from the inner shock region and diffusing into the outer 
shear layers will naturally experience an additional increase in 
energy. In 3d-hydrodynamical simulations (Aloy et al. 1999), 
the shear layer is found to broaden with distance along the jet, 
eventually becoming comparable to the beam width. If one thus 
takes Ar ~ 1 kpc near to the jet end (we note that the full width 
of the turbulent shear layer may be even larger than that seen 
in the radio), the classical Hillas confinement condition would 
constrain maximum protons energies to be E ^ 10 20 (B/10 _4 G) 
eV. This value might be somewhat higher, if the magnetic field 
is amplified in the shear (e.g., due to Kelvin-Helmholtz-driven 
shear instabilities, see Zhang et al. 2009). On the other hand, for 
smaller transverse shear widths, non-gradual shear acceleration 
(Ostrowski 1998; Rieger & Duffy 2004) is expected to take over 
earlier, leading to similar maximum energy considerations but 
different spectral characteristics. 

4. Conclusions 

If Cen A hosts a radiatively inefficient inner accretion disk, ef- 
ficient electron acceleration in its BH-jet magnetosphere and in- 
verse Compton upscattering of ambient disk photons to the TeV 
regime appears to be possible. After M87, Cen A could thus be 



another primary candidate source in which VHE processes close 
to the black hole may allow a fundamental diagnosis of its im- 
mediate environment, offering important insights into the central 
engine in AGNs. 

Whether Cen A indeed proves to be an UHECR emitting source 
may require confirmation by additional observations. From a 
theoretical point of view and based on conventional accelera- 
tion concepts, it seems challenging to account for a possible pro- 
duction of UHECR protons beyond 10 19 eV. On the other hand, 
as shown above, shear acceleration along the large-scale jet in 
Cen A may help to overcome this problem by increasing the en- 
ergy of shock-accelerated seed particles by a factor of some tens. 
If this would be the case, spectral changes in the cosmic ray spec- 
trum might be partly due to the operation of a new acceleration 
mechanism and not simply propagation effects. 

Acknowledgements. Discussions with B. Reville, A. Taylor, C. Fendt and V. 
Beskin, and useful comments by the referee are gratefully acknowledged. 
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